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ABSTRACT 

We point out a strong time-evolution of the mass-to-light conversion factor 77 commonly 
used to estimate masses of unresolved star clusters from observed cluster spectro-photometric 
measures. We present a series of gas-dynamical models coupled with the Cambridge stellar 
evolution tracks to compute line-of-sight velocity dispersions and half-light radii weighted by the 
luminosity. We explore a range of initial conditions, varying in turn the cluster mass and/or 
density, and the stellar population's IMF. We find that 77, and hence the estimated cluster mass, 
may increase by factors as large as 3 over time-scales of 50 million years. We apply these results 
to an hypothetic cluster mass distribution function (d.f.) and show that the d.f. shape may be 
strongly affected at the low-mass end by this effect. Fitting truncated isothermal (Michie-King) 
models to the projected light profile leads to over-estimates of the concentration parameter c of 
5c w 0.3 compared to the same functional fit applied to the projected mass density. 



Subject headings: stellar dynamics, stars: evolution 

1. Introduction 

The formation of star clusters in bursts of star 
formation during galactic mergers has attracted 
much attention since the ground-breaking study 
by Schweizer (1986). Young clusters can account 
for ~ 20% of the UV light flux of their sample 
starburst galaxies (compared with < 1% for the 
Milky Way; e.g. Meurer et al. 1995). Such num- 
bers bring to focus the role that star formation in 
clusters plays in shaping the overall (galactic) stel- 
lar mass function. Proto-typical cases where clus- 
ter formation has been a spectacular manifestation 
of interaction-induced starbursts are the merging 
systems NGC 4038/39 (the Antenna?) and the 
nearby galaxy M82 (recent interaction with M81). 
Many of the brightest clusters in these galaxies 
have estimated ages on the order of 10 7 years, 
based on optical and near-IR spectra. High res- 
olution spectroscopic data and HST images have 
been used to measure velocity dispersions and esti- 
mate virial masses. The line-of-sight velocity dis- 
persion aid (LOSVD) relates to the mass M and 
projected half-light radius r^ p of a cluster in virial 
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where 77 is a dimensionless free parameter. A num- 
ber of authors have set i] w 10 in their studies to 
derive M from (JTJ (Sternberg 1998; Mengel et al. 
2002; Smith & Gallagher 2001; McCrady et al. 
2003; Maraston et al. 2004). McCrady et al. give 
a derivation of this value of rj. Mengel et al. (2002) 
quote a range from 5.6 — 9.7 for King (1966) mod- 
els with concentration parameter in the range 0.5 
to 2.5, which corresponds to most galactic glob- 
ular clusters. We discuss the value of 77 further 
below. Dynamical masses derived from (fTJ have 
been compared to the stellar masses of synthetic 
populations, using both standard (field) and non- 
standard stellar initial mass functions (henceforth 
IMF; see Kroupa 2002 for a review). The data 
were found to be inconsistent with a universal IMF 
(Mengel et al. 2002; Smith & Gallagher 2001). In 
particular, several clusters in M82 were found to 
be over-luminous with respect to their estimated 
mass (low mass-to- light ratio, M/L*). This sug- 
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Fig. 1. — Parameter r\ (top lines) and half-light 
radius rh p (bottom lines with symbols) versus 
time. The initial models were unsegregated King 
^ jo 1 = 6 models. Results are shown for mod- 
els with different initial surface density (cf. text). 
When the mass range is increased from 17 M 
to 50 M Q , r\ rises sharply over the first 10 Myr 
(open circles, case 'M82-F'); thereafter it rejoins 
the curve displayed. 



gests that M82 clusters may form with a top-heavy 
stellar IMF (Smith & Gallagher 2001; McCrady et 
al. 2003). 

2. Setup 

Much of what we report in this section draws 
from Boily et al. (2005). 

2.1. Numerical method 

The equations of motion were integrated nu- 
merically based on the gas-dynamical approach pi- 
oneered by Larson (1970) and developed further 
by Louis & Spurzem (1991) to include anisotropic 
velocity fields. The method leans on an analogy 
between exchange of kinetic energy through star- 
star interactions and the classical heat-diffusion 
process of fluid dynamics (Lynden-Bell & Eggle- 



ton 1980; Heggie & Ramamani 1988). The im- 
plementation in spherical coordinates we used is 
largely due to Louis & Spurzem (1991) and Giersz 
& Spurzem (1994) but extended to include a spec- 
trum of stellar masses (Spurzem & Takahashi 
1995). 

The mass spectrum is sampled at constant log- 
arithmic increments in the interval {too, mi}; we 
used 14 bins in our standard runs ((51nm/M « 
0.329). To each mass bin corresponds the same 
initial (continuous) radial density profile. Star- 
star interactions (including those between stars of 
the same mass bin) lead to the diffusion of kinetic 
energy. Roughly speaking, the resulting change in 
the velocity dispersion of each mass bin causes a 
readjustment of the density profile. This is ob- 
tained using a semi-implicit Henyey integration. 
The gravitational potential is then updated by ap- 
plying Poisson's equation (see Louis & Spurzem 
1991). 

We checked that the mass-segregation time ob- 
tained from the models is in quantitative agree- 
ment with N-body and Fokker-Planck integrators 
by evolving Plummer models with three species 
of stars and following constant-mass Lagrangian 
radii in time for all three components (cf. Spitzcr 
& Shull 1975, Fig. 1). Spurzem & Takahashi 
(1995) report excellent agreement from their own 
two-component test calculations. 

2.2. Stellar evolution 

The different mass components are evolved 
according to the Cambridge tracks (Pols et al. 
1998; Hurley et al. 2000). The tracks are effi- 
ciently coded in the form of fitting functions of the 
kind first presented by Eggleton, Tout & Fitchctt 
(1989). The functions return the current bolomet- 
ric luminosity, radius, mass and metal abundance 
for given time and initial metal abundance Zq; we 
set Z = 0.02 (solar abundance) throughout. 

A filter can be applied to the bolometric lu- 
minosity from model stellar atmospheres, and the 
total flux in a specified waveband read by interpo- 
lation polynomials from the Basel stellar library 
(Lejeune et al. 1998). We have mapped the stellar 
luminosity near the strongest near-IR CO band- 
head (A ~ 2.2 - 2.29 fxm) and the Call triplet 
(A ~ 8200 - 8600 A) together with the bolometric 
limit (all wavelengths). Since our results are not 
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strongly dependent on waveband issues, we will 
discuss only quantities computed from bolometric 
fluxes. 

2.3. Models 

All models were scaled to N = 500, 000 mem- 
ber stars and a total mass of M ~ 2 x 10 5 M Q 
when an Salpeter mass function is implemented. 
This is significantly less than what is measured 
for massive clusters in starburst galaxies, where 
membership may well exceed 10 6 . This point is 
important because the time-scale for segregation 
scales with AT/ In AT. However we may explore 
the importance of varying the dynamical evolu- 
tion time by constructing models with different 
sizes and central densities, while keeping the to- 
tal mass, and N, constant. This approach is at- 
tractive because two systems with the same mean 
density will have the same dynamical time t CT , in- 
dependently of their mass and size. We first setup 
three models with identical if? /a 2 — 6 King param- 
eter but each with different half-light radius. The 
model labeled 'M82-F' has a mean surface density 
ks 1.5 X 10 4 M Q /pc 2 or half the value we derive for 
that cluster from the data of Smith & Gallagher 
(2001). The densest model is labeled 'R136', in 
reference to the 30 Doradus cluster (central vol- 
ume density ~ 10 6 M Q /pc 3 ; Brandl et al. 1996). 
Finally, a third low-density model was evolved for 
comparison (labeled 'Low'). 

3. Results & Applications 
3.1. Reference models 

We start with the densest model 'R136' for 
which we compute the shortest relaxation time 
i r h- Since M, rh p and (u\) are all known from 
the simulations, we solve for r\ directly from |T]). 
The run of r\ in time is displayed on Fig. [TJ The 
rapid increase of 77 from an initial value ks 8.2 is 
striking. After 15 Myr of evolution, rj has more 
than doubled. The projected half-light radius rh p 
is displayed alongside rj. The radius r\ ip decreases 
steadily in time, a direct result of the migration of 
massive stars toward the centre. Note the slow but 
systematic rise of rj after w 35 Myr, when r\ > 20 
(factor > 2.5 from its initial value. The results are 
robust and can be extended to other stellar IMF 
(Fleck et al. 2006). 



The time-evolution of rj for all three models 
is displayed on Fig. 1. Both 77 and r^ p for the 
'Low' model stay essentially constant through- 
out. However the case 'M82-F' (initial surface 
density > 10 M Q /pc 2 ) shows unmistakable evolu- 
tion, suggesting that for clusters of such or higher 
mean surface density we may no longer presume a 
time-independent mass-to-light parameter rj. This 
has to be interpreted with care since unresolved, 
bright clusters will be more massive than the ones 
modelled here; despite this caveat, evolution will 
always take place in the central regions where the 
density peaks. This has immediate bearing on fit- 
ting cluster light profiles. 

The LOSVD changes relatively little over time 
in comparison: we measure a monotonic decrease 
of (a x ) from w 10.4 to 9.3 km/s (or, -10.6%) for 
the system as a whole, although for individual 
components evolution was more significant: down 
« 30% for the most massive stars, while the light- 
est component enjoys an increase of a comparable 
magnitude. These effects can all be traced back 
to the dynamical mass segregation. 

3.2. Embedding in a merger simulation 

In this section we ask what effect the back- 
ground galactic potential may have on the forma- 
tion and early evolution of a star cluster. More 
specifically, we ask whether the tidal field is an 
ingredient that matters when it comes to under- 
standing cluster formation in mergers. The results 
reported here are taken from a doctoral disserta- 
tion by Fleck (2007). 

We built two equal mass bulge/disk/halo 
galaxy models to simulate an Antennae-like major 
merger with initial conditions similar to Barnes's 
(1988). We set up self-gravitating models using 
the magalie software included with the Nemqj 
stellar dynamical package (version 3.2.4; Teuben 
1995). The equations of motion were integrated 
with gyrfalcON (Dehnen 2002) using a Plummer 
softening kernel with e = 0.01 in units where the 
gravitational constant G = 1 and a numerical unit 
of time translates to 25 Myr. The total mass of 
each galaxy is equal to 7 mass units, where the ex- 
ponential disk scale- length and mass r& — Md = 1. 
The bulge is an Hernquist (1990) model with a 
scale parameter a — 1 and total mass Mb = 1. 



1 http: / /bima. astro, umd.edu/nemo/ 
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Thus our progenitor galaxies are "SO" galaxies 
with relatively warm Toomre Q = 1.5 disk pa- 
rameter to ensure local stability to fragmentation 
modes. Note that our model galaxies differ mostly 
from Barnes's in that he used a lighter King (1966) 
model for the bulge. The dark matter halo was 
generated from a truncated isothermal distribu- 
tion as described in Hernquist (1993). The mass 
ratio between individual particles ~ 1 to minimize 
two-body heating; a total of 700, 000 particles 
were used for each galaxy, for a mass resolution 
of (1.15 - 2.85) x 10 5 M for the different com- 
ponents. Converted to physical units, the simula- 
tions matched well the kinematics of the Antennae 
(peak systemic velocity of 130 km/s) but for a to- 
tal baryonic mass of ~ 9.2 x 10 10 M©, some twice 
that measured for that system (cf. Hibbard et al. 
2001). 

The galaxies were set on a prograde encounter 
with tidal tails developing on cither side of 
the orbital plane, implying non-coplanar disks. 
Fig. V.17 of Fleck (2007, p.154) shows the orbit 
of the centers of mass of the two galaxies together 
with the orbit of the associated point-mass prob- 
lem initially on a bound but eccentric orbit. The 
extent and mass of the dark matter halos are such 
that the galaxies are braking after the first passage 
(t ks 11) before falling back together and merging. 
Tidal tails form soon after the first passage and 
continue to expand until the end of the simulation 
at t ks 30 (for a graphical time-sequence, see Fleck 
2007, Fig. IV.9, p.123). 

3.3. Individual orbits 

Our mass resolution of » 2 x 10 5 M Q implies 
that individual disk and bulge particles represent 
cluster-size star forming regions, which we follow 
on their orbit. To explore a full range of formation 
histories, we chose several orbits that ended up ei- 
ther in the tails or remained close of the nuclear re- 
gion at the end of the simulation. We retrieved the 
tidal tensor at the particle position using a second 
order finite difference scheme from the forces at six 
points distributed on a box of radius ~ e (e = 40 
pc) . The forces were directly summed over all the 
particles. We then computed the eigenvalues {A,} 
of the tidal tensor using the fortran routines of 
Kopp (2008) and looked for regions that undergo 
purely compressive modes (max{Ai} < 0). On fig- 
ure[5J we graph a selection of orbits and their max- 



imum eigenvalue as a function of time. Note that 
the orbit going back into the nuclei (dot-dashed 
line) experiences the strongest tidal force espe- 
cially near the second passage, at and around time 
t w 23. Nevertheless, this orbit also experienced 
a fully compressive mode for a cumulative time 
interval of one unit. The orbit which spends the 
longest total period in a compressive mode (St ~ 
8 time units, or ~ 100 Myr, solid line) is the one 
staying closest to the nuclei after a short trek in 
the tidal arms. The other two orbits also expe- 
rience compressive modes immediately after the 
first passage, between t = 12 and t = 18. 

Fig. 1 of Mengel et al. (2005) shows that young 
massive clusters are dispatched throughout the re- 
gion surrounding the tidal bridge linking the two 
nuclei. To determine whether the distribution in 
space of young clusters may coincide with the com- 
pressive modes in the current-day configuration, 
we isolated a series of bodies that ended in the 
overlap region of the merger simulation and found 
that in several instances the bodies were still em- 
bedded in a volume of compressive tidal forces. Al- 
together we find that the fraction of bodies expe- 
riencing a compressive tidal field during the simu- 
lation rises from ~ 2% in the initial configurations 
at rest, to ~ 15% during both the first and second 
passage. It is clear from the statistics that tidal 
fields help push in and hold together the material 
that soon will burst to form a large ensemble of 
clusters (Renaud et al. 2008). This gathering of 
mass has long been evoked as a possible mecha- 
nism of dwarf galaxy formation with gas dissipa- 
tion (e.g., Bournaud et al. 2006). Here we see that 
the same mechanism is operative at the very heart 
of a galaxy merger on scales that match those of 
rich star formation regions. 

CMB is grateful to the organizers for an invita- 
tion to present the results of this study. 
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Fig. 2. — a) Left-hand panel: Four orbits up to the final simulation time. All particles inside a radius r < 3 
have been erased to avoid overcrowding, b) Right-hand panels: Time-evolution of the maximum eigenvalue 
measured along each orbit. Note the logarithmic scale. The top panel shows positive values, the lower panel 
negative values of compressive tides. 
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